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In Amin and Choi [l|, we show that an adiabatic quantum algorithm for the NP-hard 
maximum independent set (MIS) problem on a set of special family of graphs in which 
there are exponentially many local maxima would have the exponentially small minimum 
spectral gap and thus would require the exponential time, due to the first order quantum 
phase transition (FQPT). The problem Hamiltonian of the adiabatic quantum algorithm 
for MIS is based on the reduction to the Ising problem and has flexible parameters. In 
this paper, we show numerically on the 15-vertex graph that by choosing the parameters 
appropriately in the problem Hamiltonian (without changing the problem to be solved) 
for MIS, we can prevent the FQPT and drastically increase the minimum spectral gap. 
The result is further supported by visualization from the Decomposed State Evolution 
Visualization (DeSEV) — a visualization tool we introduced. Furthermore, our result 
also serves to concretely clarify that it is not sufficient to consider one specific problem 
Hamiltonian for proving the failure of adiabatic quantum optimization for a problem, as 
explained in [2] . We also raise the basic question about what the appropriate formulation 
of adiabatic running time should be. 

1 Introduction 

Adiabatic quantum computation (AQC) was proposed by Farhi et al. 0] in 2000 as an 
alternative quantum paradigm to solve NP-hard optimization problems, which are believed 
to be classically intractable. The same idea to the adiabatic quantum optimization, under a 
different name of quantum annealing, was first put forward by Apolloni et al. in 1988, see 016] 
and references therein for a history of the field. It was shown that AQC is not just limited to 
optimization problems, and is polynomially equivalent to conventional quantum computation 
(quantum circuit model) [3 [8]. A quantum computer promises extraordinary power over 
a classical computer, as demonstrated by Shor [9 in 1994 with the polynomial quantum 
algorithm for solving the factoring problem, for which the best known classical algorithms are 
exponential. Just how much more powerful are quantum computers? In particular, we are 
interested in whether an adiabatic quantum computer can solve NP-complete problems more 
efficiently than a classical computer. 

Unlike classical computation or quantum circuit model in which an algorithm is specified 
by a finite sequence of discrete operations via classical/quantum gates, the adiabatic quan- 
tum algorithm is continuous. It has been assumed (see Section^ for more discussion) that, 



according to the adiabatic theorem, the dominant factor of the adiabatic running time (ART) 
of the algorithm scales polynomially with the inverse of the minimum spectral gap g m \ n of 
the system Hamiltonian (that describes the algorithm). Therefore, in order to analyze the 
running time of an adiabatic algorithm, it is necessary to be able to bound g mm analytically. 
However, <7 m j n is in general difficult to compute (it is as hard as solving the original problem if 
computed directly). Rigorous analytical analysis of adiabatic algorithms remains challenging. 
Most studies have to resort to numerical calculations. These include numerical integration 
of Schrodinger equation [H [TO] , eigenvalue computation (or exact diagonization) [TTJ [T2] , and 
quantum Monte Carlo (QMC) technique [T31 H3] - However, not only are these methods lim- 
ited to small sizes (as the simulations of quantum systems grow exponentially with the system 
size), but also little insight can be gained from these numbers to design and analyze the time 
complexity of the algorithm. 

Perhaps, from the algorithmic design point of view, it is more important to unveil the 
quantum evolution black-box and thus enable us to obtain insight for designing efficient adia- 
batic quantum algorithms. For this purpose, we devise a visualization tool, called Decomposed 
State Evolution Visualization (DeSEV). 

One of the original ideas of AQC in [3113] was proposed as an energy minimization algorithm 
that aims to use coherent quantum evolution to avoid trapping in the local minima that trip 
classical algorithms of NP-hard optimization problems. However, several works [T3 E3 IT7] 
showed that their adiabatic quantum algorithm failed to avoid local minima and would take 
exponential time for some problem. As we pointed out in [2], their lower bounds are only for 
their specific adiabatic quantum algorithms, and their arguments are not sufficient for showing 
the lower bound of all adiabatic quantum algorithms of the same problem. Nevertheless, one 
still might argue that their results provide "convincing evidence" that AQC would fail to solve 
problems with many local minima. For this purpose, we set out to study graphs that have 
exponential many local maxima for the NP-hard Maximum Independent Set (MIS) problem 
which AQC "naturally" solves, see Section [3] and Choi [18]. While it is not difficult to come 
up with such graphs, the challenge lies in coming up such graphs with small sizes such that we 
can visualize the quantum evolution using DeSEV, which relies on the (numerical) eigenvalue 
computation. After many trial-and-errors, with the aid of DeSEV, we constructed a special 
family of graphs in which there are exponentially many local maxima for MIS, with the 
smallest size of such graph being 15. For the reference sake, we call this family of graphs 
the CK graph. The numerical results of an adiabatic algorithm for MIS on these graphs 
suggested that g mm is exponentially small and thus the algorithm requires exponential time. 
These results were then explained by the first order quantum phase transition (FQPT) by 
Amin and Choi in [JJ. That is, our result agreed with the speculation — the system got 
trapped in the local minima and the particular AQO algorithm failed due to the FQPT. 
Since then, there have been some other papers (Altshuler et al., [30] ; Farhi et al., |21) ; 
Young et al., 14J; Jorg et al., [35J[53]) investigating the same phenomenon, i.e., first order 
quantum phase transition. In particular, Farhi et al. in [2T] suggested that the exponential 
small gap caused by the FQPT could be overcome (for the set of instances they consider) 
by randomizing the choice of initial Hamiltonian. In this paper, we show numerically that 
by changing the parameters in the problem Hamiltonian (without changing the problem to 
be solved) of the adiabatic algorithm for MIS on the CK graph, we prevent the FQPT from 
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occurring and significantly increase g m in . We further support our result by visualization from 
DeSEV. To the best of our knowledge, this is the first time that such a numerical result 
has been shown. We do so by scaling the vertex-weight of the graph, namely, multiplying 
the weights of vertices by a scaling factor. In order to determine the best scaling factor, we 
raise the basic question about what the appropriate formulation of adiabatic running time 
should be. We remark that the scaling factor in turn relates to the bit of precision required 
for the parameters in the Hamiltonian which is an important physical resource and plays a 
critical role in the computational complexity. Finally, our result serves to further clarify (see 
[H I19j for explanation) that it is not sufficient to consider one specific problem Hamiltonian 
for proving the the failure of adiabatic quantum optimization for a problem. 

This paper is organized as follows. In Section [2j we review the adiabatic quantum algo- 
rithm, and the adiabatic running time(ART). In Section |3j we recall the adiabatic quantum 
algorithm for MIS based on the reduction to the Ising problem. In Section [4j we describe 
the visualization tool DeSEV and the CK graph. We show examples of DeSEV on the MIS 
adiabatic algorithm for the CK graph. In Section [5j we describe how changing the parameters 
affects <? m in, and raise the question about ART. We conclude with the discussion in Section [6] 

2 Adiabatic Quantum Algorithm 

An adiabatic quantum algorithm is described by a time-dependent system Hamiltonian 

H(t) - (1 - s(t))H- tnit + s(i)H prob lem (1) 

for t £ [0, T], s(0) = 0, s(T) = 1. There are three components of %{■): (1) initial Hamil- 
tonian: H(0) — H\ n \t', (2) problem Hamiltonian: H(T) — ^problem! and (3) evolution path: 
s : [0,T] — > [0, 1], e.g., s(t) = |r. T~L{t) is an adiabatic algorithm for an optimization problem 
if we encode the problem into the problem Hamiltonian "Hprobiem such that the ground state 
of ^problem corresponds to the answer to the problem. 

In this paper, we fix the evolution path by the linear interpolation function s(t) = Jp. 
Hereafter, we describe an adiabatic algorithm by the re-parametrized Hamiltonian 

U( s ) = (1 - s )^init + sHproblem (2) 

where s £ [0,1], with understanding that s(t) = t/T. For a more general interpolation 
path see 24 . Furthermore, throughout this paper, we fix the initial Hamiltonian to be 
H\ nlt = — Sjgv(G) a i- When it is clear from context, we also refer to the problem Hamiltonian 
as the adiabatic algorithm for the problem. 

2.1 Adiabatic Running Time 

In their original work [3] , the running time of the adiabatic algorithm is defined to be the same 
as the adiabatic evolution time T, which is given by the adiabatic condition of the adiabatic 
theorem. However, this definition is under the assumption of some physical limit of the 
maximum energy of the system (see e.g., |25j). and is not well-defined from the computational 
point of view, as observed by Aharonov et al. [7]. They re-define ART('H) as T-max s ||%(s)||, 
taking into the account of the time-energy trade-off in the Schrodinger's equation^ 1 

a Namely, i d ^ ]) = T ■ H(s)\ip(s)) = § • KH{s)\ip(s)) where K > is a constant. 
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On the other hand, given the extensive work on the rigorous proofs of the adiabatic 
theorem, it is interesting (if not confusing) that many different versions of the adiabatic 
conditions have been recently proposed. These include ^71 1^1 15T | [5^ 1 155 1 [M l [55 1 155 1 157] 
in the quantum physics community, and (T7J [7J in the computer science community. Most 
of these studies imply that ART scales polynomially with the inverse of the spectral gap of the 
system Hamiltonian, which is sufficient when one is interested in the coarse computational 
complexity of algorithms, namely, the distinction between polynomial and exponential running 
time. 

However, from both the practical and algorithmic point of view, it is important to have 
a more precise formulation of ART. First, this is because the specification of the adiabatic 
evolution time T is required in an adiabatic algorithm, and therefore a tight and simple upper 
bound is desired. That is, while there are complicated formulas such as the ones from [37], 
although accurate, they are not useful if the formulas can not be efficiently evaluated. Second, 
we are interested in the actual time complexity of the algorithm, and not just the polynomial 
vs. exponential distinction. It is necessary to have a more precise formulation of ART such 
that basic algorithmic analysis can be carried out. Third, at this stage of research, it is 
particularly important to have such a formulation because the spectral gap, which plays the 
dominating role in the formulation of ART, is difficult to analyze. All current efforts on the 
spectral gap analysis resort to numerical studies, and that means the studies are restricted 
to small problem sizes only. Therefore, to gain insight into the time complexity of algorithms 
from these small instances, it is important that the formulation of ART applies to small sizes. 
So what is the appropriate formulation of ART? What should the adiabatic condition(s) be? 
This is in contrast to the study in [23] where the exact form of adiabatic condition is not 
essential. In Section |5.2[ we compare three closely related versions and raise the question 
about what the appropriate adiabatic running time should be. 

3 An Adiabatic Algorithm for MIS 

In this section, we recall the adiabatic algorithm for MIS that is based on the reduction 
to the Ising problem, as described in [TH]. First, we formally define the Maximum- Weight 
Independent Set (MIS) problem (optimization version): 

Input: An undirected graph G — (V(G),E(G)), where each vertex i £ V(G) = 
{1, . . . , n} is weighted by a positive rational number Cj 

Output: A subset S C V(G) such that S is independent (i.e., for each i,j £ V(G), 
i 7^ i) ij E(G)) and the total weight of S ~ J^ies c i lB maximized. Denote the optimal set 
by mis(G). 

There is a one-one correspondence between the MIS problem and the Ising problem, which 
is the problem directly solved by the quantum processor that implements 1/2-spin Ising 
Hamiltonian. We recall the quadratic binary optimization formulation of the problem. More 
details can be found in [15] , 

Theorem 1 (Theorem 5.1 in |18j ) // Jij > min{ Ci,Cj) for all ij £ E(G), then the 
maximum value of 




(3) 



iGV(G) 



ijeE(<?) 
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is the total weight of the MIS. In particular if Jij > min{cj,Cj} for all ij £ E(G), then 
mis(G) = {i £ V(G) : x* = 1}, where (x*, ...,<) = argmax {xi ^,^ Xn)e{0A}n y(x 1 , . . .,x n ). 

Here the function y is called the pseudo-boolean function for MIS. Notice that in this 
formulation, we only require > minjcj, Cj}, and thus there is freedom in choosing this 
parameter. In this paper we will show how to take advantage of this. 

By changing the variables (xt — ^r 1 ), it is easy to show that MIS is equivalent to 
minimizing the following function, known as the Ising energy function: 

£(si,...,s n ) = ^2 h t Si+ ^2 JijSiSj, ( 4 ) 

»eV(G) ijeE(G) 

which is the eigenfunction of the following Ising Hamiltonian: 

"Hlsing = ^2 hiCT i + X! Jij^iVj ( 5 ) 
»GV(G) »jeE(G) 

where hi = Ejenbr(i) J u" _ 2c ^ nbr (*) = : *.? € E (G)}, for i £ V(G). 
4 DeSEV and CK Graph 

In this section, we describe a visualization tool, called Decomposed State Evolution Visual- 
ization (DeSEV), which aims to "open up" the quantum evolution black-box from a compu- 
tational point of view. Consider the above adiabatic algorithm for MIS. Recall that according 
to the adiabatic theorem, if the evolution is slow enough, the system remains in the instan- 
taneous ground state. Let \ip(s)} be the ground state of V.(s), for s £ [0, 1]. For a system of 
n-qubits, \4>(s)) is a superposition of 2" possible computational states, namely, 

l^( s )) = Oi x (s)\x), where ^ \a x (s)\ 2 = 1. 

rrG{0,l}" x£{0,l}" 

For example, we have the initial ground state \ip{0)) = J2xe{o i}™ \ x )> which is the uniform 
superposition of all 2™ states, while the final ground state IV'(l)) = |#i#2 • • ■ x Vli corresponding 
to the solution state. A natural question is: what are the instantaneous ground states \ip{s)), 
for < s < 1, like? In particular, we would like to "see" how the instantaneous ground 
state evolves? A naive solution would be to trace the 2™ amplitudes a x . The task becomes 
unmanageable even for n = 10, which has 1024 amplitudes, even though many may be 
negligible (close to zero). 

To make the "visualization" feasible, we introduce a new measure IV Suppose that H(l), 
has (m + 1) < 2™ distinct energy levels: E < E\ < ... < E m . For < k < m, let 
Dk = {x £ {0, 1}™ : H(l)\x) = Ek\x)} be the set of (degenerate) computational states that 
have the same energy level Ek (with respect to the problem Hamiltonian H(l)), and define 

r k (s)= J2 KWI 2 - 

xeD k 

In other words, Tfe(s) is the total percentage of (computational) states of the same energy 
level Ek participating in \ip(s)). The idea is now to trace Tk instead of a x . Here we remark 
that Tfe are defined for any eigenstate and not just for the ground state. 
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For our purpose, we constructed a special family of vertex-weighted graphs for the MIS 
problem. We designed the problem instances such that the global minimum is "hidden" in 
the sense that there are many local minima to mislead local search based algorithms. Note 
that the size of the smallest instances needs to be necessarily smaller than 20 as we are relying 
on the eigenvalue computation (or exact diagonization) to compute T^. 

CK Graph Construction. Let r, g be integers, and wa, wb be positive rational numbers. 
Our graphs are specified by these four parameters. There are two types of vertices in the 
graph: vertices of a 2g-independent set, denoted by Va, and vertices of g r-cliques (which 
form r 9 maximal independent sets), denoted by Vb- The weight of vertex in Va (Vb resp.) is 
wa {wb resp.). The graph is connected as follows. Partition the vertices in Va into g groups 
of 2 (independent) vertices. There are also g groups of r-cliques in Vb- We label both groups 
accordingly such that each group in Va is adjacent to all but one (the same label) r-cliques in 
Vb- Note if wb < 2wa, then we have Va forming the (global) maximum independent sets of 
weight 2gwA, while there are r 9 (local) maximal independent set of weight gins- See Figure 
[T] for an example of a graph for r = 3 and g = 3. In general, there are infinitely many such 
graphs specified by the parameters r and g. 




(a) (6) 

Fig. 1. (a) A CK graph for r = 3 and g = 3. The graph consists of 15 vertices: Va = {1, ■ ■ • , 6} 
forms an independent set of size 6, while Vb, consisting of g(= 3) groups of r(= 3) triangles: 
{7, 8, 9}, {10, 11, 12}, and {13, 14, 15}, forms 3 3 independent sets of size 3. The graph is connected 
as follows. The 6 vertices in Va are divided into 3 groups: {1, 2}, {3, 4}, and {5, 6}. The vertices in 
each group are adjacent to vertices in two groups of three triangles in Vb (as illustrated by different 
colors), (b) The drawing of the graph with explicit connections. The weight of a vertex in Va 
(Vb resp.) is wa (wb resp.). We set wa = 1, and consider 1 < wb < 2. For explanation purpose, 
we represent a vertex in Va by a •, and a vertex in Vb by a A. Therefore, Va = {•, •, «, •, }, 
forms the MIS of weight 6; while {A, A, A} is a maximal independent set of weight 3wb(< 6). 



4-1 DeSEV for the MIS Adiabatic Algorithm on a 15-vertex CK Graph 

In the section, we fix the CK graph with r = 3, g = 3 as illustrated in Figure [l] We set 
u>a = 1, and consider 1 < wb < 2. The graph G consists of 15 vertices: Va = {1, ... ,6} 
forms the maximum- weight independent set of weight 6; while Vb, consisting of 3 groups of 3 
triangles: {7, 8, 9}, {10, 11, 12}, and {13, 14, 15}, forms 3 3 maximal independent sets of weight 
3wb < 6. 

According to Eq.Q, the problem Hamiltonian (and thus the adiabatic algorithm) for MIS 



G 



^(6J-2)a?+ ^(6J-2w B )o? + J J] 

iev A iev B yeE(G) 

2 > w B for all ij € E(G). 

Notation on representing the computational states. For a computational state |sci£2 ■ ■ ■ x n ) 
where Xi € {0, 1}, we adopt the zero position representation, namely, represent it by \i1i2 ■ ■ - ik) 
where Xj = if and only if j — i t for some t. That is, we represent 1 000000 111111111} (the 
solution state) by 1 123456). Further, we use a • to denote a vertex in Va, a A for a vertex in 
Vb- That is, the solution state is now represented by (••••••}, while |AAA), corresponding 

to a local maximal independent set of weight 3wb with one vertex from each triangle. 

Maximum vs Minimum. The maximum of MIS corresponds to the minimum of the 
Ising energy. For explanation purpose, instead of referring to the energy values of the Ising 
Hamiltonian, we will refer to the values of MIS given by the pseudo-boolean function y in 
Eq.([3| by "(-)energy", where "(-)" is to indicate the reverse ordering. 

Example. The (-)energy of !••••••) is 6; while |AAA— A) is Awb — J, where A— A 

represents two connected vertices from Vb, e.g. vertex 7 and 8 in Figure [T] 

See Figure [2] for the DeSEV of the the ground state of the adiabatic algorithm with Hi 
in Eq.([6]) as the problem Hamiltonian for wb — 1.5 and 1.8. 

4-2 FQPT and Perturbation Estimation 

To gain better understanding, in [1], we vary the weights of vertices: fix wa = 1, while varying 
u>b from 1 to 1.9 with a step size of 0.1. That is, we fix the global maximum independent 
set, while increasing the weight of the local maximum. As the weight of wg increases, the 
minimum spectral gaps get smaller and smaller (indeed, from 10 _1 to 10 -8 as wb changes 
from 1 to 1.9 as shown in Table [3]). 
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Fig. 3. The minimum spectral gap <? m j n (and position s*) changes as wg changes from 1 to 1.9, 
for the (unsealed) problem Hamiltonian Hi in Eq.|6]l. 

This was consequently explained by the FQPT in [T|. By FQPT, here we mean that 
there is a level anti-crossing between two states as illustrated in Figure [4j The minimum 
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s* = 0.3805, g n 


lin = 2.04 x 10~ 2 s* 


= 0.6276, g mm = 1.04 x 10 
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Fig. 2. DeSEV (only the 7 lowest energy levels shown) of the ground state of the MIS adiabatic 
algorithm with "Hi in Eq[6]as the problem Hamiltonian for wg = 1.5 (left) and wg = 1.8 (right). 
The x-axis is the time s. The y-axis is the (-)energy level. Each color corresponds to an energy 
level. The correspondence between (-)energy levels and the states are shown. The z-axis is T. 
s* is the position of the minimum spectral gap(g m | n ). As time s increases, one can see how T of 
each energy level evolves to get some sense of the evolution. For example, for wg = 1.5 (left), for 
the (-)energy level 6 (which corresponds to the solution state), shown in brown, T changes from 
almost at s = 0.2, to more than 0.4 at s = 0.4, to almost 1.0 at s = 0.8. For wg = 1.8 (right), Y 
of (-) energy level 6 changes from almost before s = 0.6 to more than 0.9 at s = 0.7; while T of 
(-) energy level 5.4, which corresponds to the local minima, gradually increases from s = to 0.6, 
but almost after s = 0.6. 



spectral gap (g m in) and the position (s*) were then estimated based on the assumption of 
the level anti-crossing between the global minimum and the local minima using perturbation 
method. In particular, g mm was estimated by the tunneling amplitude between the global 
minimum and the local minima. The formula so derived involves combinatorial enumeration 
of the all possible paths between local minima and the global minimum, and suggested g m - ]n is 
exponentially (in terms of the problem size) small. See also [20l [1] |21j [14] for more explanation 
on the FQPT and the level anti-crossing. 

5 Varying Parameters in the Problem Hamiltonian for MIS 

In this section, we show that by changing the parameters in the problem Hamiltonian for MIS 
on the CK graph, the FQPT no longer occurs and we can significantly increase g mm . 

Recall that in the pseudo-boolean formulation of MIS as in Theorem [l] the requirement 
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(Zoom:s = 0.627 ... 0.628) 
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(-)energy 
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s* = 0.6276, g min = 1.04 x 10~ 5 

(a) (b) (c) 

3.8 4 4.8 5 5.2 5.4 6 

| • «A) • • • •> A-AA-AAA) + IAA-AA-A) | AAAA) AAA) ...... 

Fig. 4. DeSEV of the ground state and the first excited state of the MIS adiabatic algorithm with 
Hi in[(S]as the problem Hamiltonian for wg = 1.8 (a) s = . . . 1; (b) Zoom in s = 0.627 . . . 0.628; 
(c) The lowest two energy levels of H(s), s = 0. . . 1. The inset illustrates a level anti-crossing 
between two states |B) and |A), or the system has a FQPT from \B) to \A) at the anti-crossing 
s*. In this example, \A) = \ > + | > and \B) = |AavAAA) + |AAAA} + |AAA). 
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for Jij is at least min{ci, Cj}, for each ij G E(G). For simplicity, we consider the simplest 
case in which = J for all ij € E(G). In other words, we have the corresponding problem 
Hamiltonian: 

ieV(G) ij'GE(G) 

where d» is the degree of vertex i G V(G). 

The natural question is: how does the ART change when we vary J? Note that it is not 
sufficient to consider only the minimum spectral gap change (as almost all the other works 
on adiabatic quantum computation did) because by increasing J, the maximum energy of 
the system Hamiltonian also increases. Instead, in order to keep the maximum energy of the 
system Hamiltonian comparable, we keep J fixed and vary Cj instead, namely multiplying all 
weights Cj by a scaling factor, say l/k, for k > 1, which does not change the original problem 
to be solved. We remark that this is equivalent to multiplying J by k, and then multiply the 
problem Hamiltonian by (l/k). 

That is, we consider the following (scaled) problem Hamiltonian 

n k = (Jdi-2 Ci /k)a!+ Y Jct M ( 8 ) 

ieV(G) ijeE(G) 

where k > 1 is the scaling factor. 

5.1 Minimum Spectral Gap g mm Without FQPT 

The DeSEVs of Hi and W 10 are shown in Figure [5] and Figure [7] The anti-crossing between 
the global minimum and the local minima (for k = 1) no longer occurs for k = 10, and 
<7 m i n increases from 1.04 x 10~ 5 to 1.45 x 10 _1 . A worthwhile observation is the change in 
the lowest few excited energy levels: for k = 1, the lowest few excited states (beyond the 
first excited state) of the problem Hamiltonian is mainly the superposition of states from Vg 
(A) (which constitutes the local minima); while these states of the scaled (k = 10) problem 
Hamiltonian is mainly the superposition of states from Va(*) (which constitutes the global 
minimum). That is, the minimum spectral gap can be increased drastically (by as much as 
four order of magnitude in this example) when the second or higher excited energy levels are 
changed (while the lowest and first excited energy level stay the same). The DeSEVs of Hk 
for k = 1, 2, 3, 5, 10, 50 are shown in Figure [HJ 

In PP, based on the FQPT assumption, we estimate (? m j n (for "Hi) by the tunneling am- 
plitude between the local minima and the global minimum, which suggests that g mm is expo- 
nentially small. However, for k = 10, from our numerical data and DeSEV in Figure [5] we 
see that the FQPT (that causes <7 m i n to be exponentially small) no longer occurs, and g m i n 
increases significantly. This seems to suggest that g m \„ to be polynomially small instead for a 
general CK graph of size n. The problem for analytically estimating <7 m i n of T~Lk for a general 
CK graph of size n remains open. We remark here that the perturbation method is still 
valid (in fact, as we increase k, we also increase the minimum spectral gap position s* — >• 1), 
however we can no longer assume that g m \ n can be approximated by the tunneling amplitude 
between the two (localized) states. 
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5.2 Scaling Factor and ART 

In this section, we discuss what the good scaling factor should be, and how it affects the 
ART. To address this question, we need an appropriate formulation for ART. Notice that 
even for numerical studies, it is not sufficient to just consider g min (as the other numerical 
works on adiabatic quantum computation did, see e.g. |14|). but the matrix element of the 
time derivative of the Hamiltonian also matters. In particular, we adopt the following three 
formulations, which are related to the widely used traditional condition: 

f ART! (ft) = m&X °- s £ M{s) maxo^ ||ft(s)|| 
(*)< ART 2 (ft) = ^pmax < s <i||ft( S )||, where g mm = - E (s*) 

[ ART 3 (ft) = max < s <i (i ; 1 ( s ^ ) ( s)) 2 max <s<i ||ft(s)|| 

where M(s) — \(E\(s)\^\Eo(s))\ is the matrix element of the time derivative Hamiltonian 
at time s, and H(s)\Ei(s)) = Ei(s)\Ei(s)) . See Table[6]for the numerical comparisons. 
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0.54578285 


6.37e-03 


2.04c+00 


2.04e+00 


2.48e+02 


1.24c+07 


1.24c+07 


3 


0.54467568 


3.30e-02 


1.41c+00 


1.41c+00 


2.55e+02 


3.32e+05 


3.32e+05 


4 


0.55610853 


6.83e-02 


1.18c+00 


1.18c+00 


2.59c+02 


6.57e+04 


6.58e+04 


5 


0.57419149 


9.67e-02 


1.06c+00 


1.07c+00 


2.61e+02 


2.96e+04 


2.99c+04 


10 


0.66773072 


1.45e-01 


7.48e-01 


7.92e-01 


2.66e+02 


9.45e+03 


1.00e+04 


20 


0.80170240 


1.30e-01 


4.72e-01 


5.68e-01 


2.68e+02 


7.48e+03 


9.01e+03 


30 


0.99318624 


7.97e-02 


8.95e-09 


4.26e-01 


2.69e+02 


3.78e-04 


1.80c+04 


40 


0.99642154 


5.99e-02 


4.90e-10 


4.35e-01 


2.69e+02 


3.67e-05 


3.26e+04 


50 


0.99779592 


4.79e-02 


5.30e-ll 


4.41e-01 


2.69e+02 


6.20e-06 


5.16e+04 



k 


s' 


9(s') 


Mis') 


M(s') 

a(s') 2 


max s ||ft || 


ART 3 


1 


0.62763727 


1.04e-05 


4.02e+00 


3.70c+10 


2.26e+02 


8.34e+12 


2 


0.54578226 


6.37e-03 


2.04c+00 


5.02c+04 


2.48e+02 


1.24e+07 


3 


0.54461081 


3.30c-02 


1.41c+00 


1.30c+03 


2.55c+02 


3.32c+05 


4 


0.55545411 


6.83e-02 


1.18c+00 


2.54c+02 


2.59e+02 


6.57e+04 


5 


0.57223394 


9.68e-02 


1.07c+00 


1.14c+02 


2.61c+02 


2.97e+04 


10 


0.65682886 


1.46e-01 


7.75e-01 


3.64e+01 


2.66e+02 


9.66e+03 


20 


0.77115481 


1.33e-01 


5.41e-01 


3.08e+01 


2.68c+02 


8.24e+03 


30 


0.83962780 


1.08e-01 


4.43e-01 


3.82e+01 


2.69e+02 


1.02e+04 


40 


0.88050519 


8.82e-02 


3.93e-01 


5.05e+01 


2.69e+02 


1.36e+04 


50 


0.90581875 


7.39e-02 


3.63e-01 


6.64c+01 


2.69e+02 


1.79c+04 



where g(s) = Si(s) - E (s), and s' = argmax < s <i-^U. 

Fig. 6. ARTi, ART 2 , ART 3 for H k in Eq.Js). Observations: (1) q m \„ increases as k increases 
from 1 to 10, but decreases from 10 to 50. (2) ARTi, ART2, and ART3 are close for k < 5. 
(3) The matrix element M(s*) at the position of minimum spectral gap is extremely small for 
k > 30. (4) For k > 10, s* (the position of the minimum spectral gap) is different from s' , where 

/ M(s) 
s = orpmoxo< s <i^jf. 

From Table [6j we see that g m i n increases as k increases from 1 to 10, however, decreases 
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from 10 to 50 (even though it is still much larger than k — 1). The latter, perhaps, can 
be explained by the following: as k increases, the difference between the low energy levels 
decreases, and becomes dominate for k > 10. We remark that the optimal value for k seems 
to depend only on the vertex weights (for which J depends on), and independ of the problem 
size. By increasing the scaling factor, we also increase the precision (or dynamic range) 
requirement for representing the parameters (hi & Jy) in the problem Hamiltonian, which is 
one of the important physical resources. 

Notice that ART2 < ARTi < ART3. The condition given in ART3 is the formula that 
one would derive from the adiabatic approximation. The condition in ARTi is the widely 
used traditional version. The condition in ART2 was mentioned in |13j . The natural ques- 
tion is: when are they asymptotically equivalent? Indeed, the three versions of ART co- 
incide for some Hamiltonians (e.g. for k = 1). However, they can be very different for 
the large k. The main reason is that the matrix element A4(s) can be extremely small at 
the minimum spectral gap position s* . For example, for k = 50, s* — > 1, Ai(s*) is ex- 
tremely small. Note one can show that AA(s) = \(Ei(s)\H\ n \t\Eo(s)\/ s for s 6 (0, 1] because 
M(s) = 1(^(^(1) -H(0)|£ ( S )| = |(£i(s) r (s) ^ (0) |£o(s)| = K£a( S )|H(0)|£ ( S )|/ S . 
Thus, for our initial Hamiltonian, A4(s) measures the overlap of the states with one single 
bit flip, and in this case it is extremely small. Observe that the position of the minimum 
spectral gap s* is not the same as the position s' where ~^~j2 is maximized. What should be 
the appropriate formulation of ART? Should it be ART3? If so, under what condition, can 
ARTi be a good approximation to ART3? and under what condition, can we assume that 
g min is the dominating factor (as have been assumed by all other works)? 

6 Discussion 

In this paper, we have shown that by changing the parameters in the problem Hamiltonian 
(without changing the problem to be solved) of the adiabatic algorithm for MIS on the CK 
graph, we prevent the FQPT, that causes the exponential small g m \m from occurring and 
significantly increase g mm - We do so by scaling the vertex- weight of the graph, namely, multi- 
plying the weights of vertices by a scaling factor. In order to determine the best scaling factor, 
we raise the basic question about what the appropriate formulation of adiabatic running time 
should be. 

In [161115] . van Dam et al. argued that adiabatic quantum optimization might be thought 
of as a kind of "quantum local search" , and in |15j , they constructed a special family of 
3SAT instances for which the clause- violation cost function based adiabatic algorithm required 
exponential time. In [5] , we point out that the exponential small gap argument does not apply 
to a different adiabatic algorithm for 3SAT. Our CK graph was designed to trap local search 
algorithms in the sense that there are many local minima to mislead the local search process. 
From DeSEV on a 15-vertex CK graph, we see that indeed this is the case for Hi and the 
adiabatic algorithm would require exponential time due to the exponential small <7 m i n caused 
by the FQPT or the level anti-crossing between the global minimum and the local minima. 
However, for Hk (say k = 10), the FQPT no longer occurs and g mm increases significantly, 
which might suggest the possibility of exponential speed-up over "Hi. It remains challenging on 
how to analytically bound <? m i n and/or ART of the adiabatic algorithm for Hk on general (CK) 
graphs. One worthwhile observation from this work is that the minimum spectral gap can be 
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increased drastically when the second or higher excited energy states are overlapping with the 
ground state in spite of the large amount of first excited states (which constitutes the local 
minima). Recall that NP-complctc problems can be polynomial reducible to each other. The 
reduction requires only the solution to be preserved, i.e. there is a polynomial time algorithm 
that maps the solution to the reduced problem to the solution to the original problem and 
vice versa. In other words, the reduction might only preserves the solution (i.e. the ground 
state) and alter the energy levels of the problem Hamiltonian. Therefore, according to the 
observation, different reduction is possible to give rise to different problem Hamiltonians, and 
thus different adiabatic quantum algorithms, for the same problem. Indeed, we have shown 
in [5] that based on the NP-complete reductions, we describe different adiabatic quantum 
algorithms to which the arguments in |15l I20| for the failure of their adiabatic quantum 
algorithms do not apply. 

In summary, although our result is only numerical and supported by visualization, this 
small example, nevertheless, serves to show that it is possible to avoid FQPT, and also 
to concretely clarify that it is not sufficient to consider one specific problem Hamiltonian 
(and thus one specific adiabatic quantum optimization algorithm) for proving the failure of 
adiabatic quantum optimization for a problem. 
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k = 1 



k = 10 



First Excited State 



First Excited State 
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time s 
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Ground State 




time s 
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1.04 x 10" 
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AA AA) 
5/fc 
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AAA) 



6/fc 



10 3/fc 3. 6/fc 3.8/fc 4/fc 5/fc 5.4/fc 

!•••> |AA> l««A) . ... I > AAA) | ) 

Fig. 5. DeSEV of the ground state and the first excited state of the MIS adiabatic algorithm 
with problem Hamiltonian Hi (left) and Hio (right) where wb = 1-8. Notice the differences in 
the lowest few excited states. For k = 1, the 2nd and 3rd excited states are superpositions of As 
(vertex in Vb which constitutes the local optima); while for k = 10, the 2nd and 3rd excited states 
are superpositions of »s (vertex in Va which constitutes the global optimum). As a result, the first 
order phase transition from local minima to global minimum occurs for k = 1, which results in the 
9min = 1-04 X 10 -5 at s* = 0.627. For k = 10, such crossing no longer occurs, and g mln = 0.145 at 
s* = 0.667. See Figure [7] for the : 



zoom-m. 



16 



k = l(s: 0.627 . . . 0.628) k = l(s : 0.62763 . . . 0.62764) k = 10(s : 0.667 . . . 0.668) 





3.8 4 4.8 5 5.2 5.4 6 

| • «A> . ... | •••••> + |AA-AA-A) A A A-A ) AAA) ...... 

3/fc 3.6/fc 3.8/fc 4/fc 5/fc 5.4/fc 6/fc 

!•••> |AA> | ••A) . ... I > AAA) | ) 

Fig. 7. Zoom around the position s* of the minimum spectral gap. 
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k = 1 



k = 2 



k = 3 




s* = 0.5741, flf min = 9.67 x 10~ 2 s* = 0.6677, g mm = 1.45 x 10" 1 s* = 0.9977, g mm = 4.79 x 10~ 2 

k — 1 3.8 4 4.8 5 5.2 5.4 6 

| • «A> |*«««) A-A A-A ) | • • • ••> + | AA-AA-A) A A A-A ) AAA) ......} 

k as 2 1.6 1.8 1.9 2 2.5 2.7 3 

AAA) A A) • «A) !••••) •••••) AAA) ...... 

k > 3 3/fc 3.6/fc 3.8/fc 4/fc 5/fc 5.4/fc 6/fc 

!•••> |AA> | ••A) ••••> I > AAA) | > 

Fig. 8. DeSEV of the ground state and the first excited state of the adiabatic algorithm with 
problem Hamiltonian Hk for wb = 1.8, where k = 1,2,3,5,10,50. s* is the position of the 
minimum spectral gap(g m ; n ). 
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